PCA analysis of ITD auditory evoked field 
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1 Introduction 

Distribution of MEG (Magnetoencephalography) 
data varies due to the number of the event factors in 
stimuli, due to the contributions of neuronal 
activities evoked by those event factors. 
Characterizing the distribution in spatial and 
temporal ways would be helpful for analyzing 
biomagnetic data containing the contributions of the 
activities related to those event factors. In this paper, 
we describe the spatial and temporal 
characterizations of the distribution using PCA 
(Principal Component Analysis). First, the spatial 
distribution of the data was characterized by 
comparing the curve of the cumulative contribution 
ratio. Second, temporal distribution is characterized 
by finding contribution of the specific biomagnetic 
activity in the data using an extension of MUSIC 
(Multiple Signal Classification) [1]. 

2 Methods 

2.1 Spatial characterization 

PCA is employed to spatially characterize the 
distribution of the data. MEG data is expressed by 
x(^) , and number of sensor channels is M. 
Decomposition of the covariance matrix of x(f) is 
expressed as, 

(x(0x( 0 T ) = VAV t , (1) 

where denotes time average, 

A = diag(A^, ? (A^ — A 2 — ••• — A^) ? and 

V = (v 1? v 2 ,...,v M ) , ( V, is the eigenvectors 
corresponding to A.). The cumulative contribution 
c(K ) ratio up to Kth component is expressed as, 

c(K) = £*,/£ A (2) 

z=l / 7=1 

2.2 Temporal characterization 

To find the contribution of a specific biomagnetic 
activity in MEG data, we propose the virtual 
beamformer method. MEG data having spatial and 
temporal pattern of the activity in interest is selected 


as the specific activity data a (t) . By applying PCA, 
the matrix of eigenvalues A and the matrix of 
eigenvectors V are obtained as shown in Eq. 1. 
Considering the number of components K that c(K ) 
in Eq. 2 is no longer increasing, the eigenvectors in 
V can be classified with signal and noise 
components so that the corresponding subspaces are 
spans of V signal = (v,, v 2 ,..., v *) and 

Vioise — ( V AT+1 > V A:+2 ***** ) * 

Now, if b(0 is the target data to evaluate, then the 
contribution score e{t) of a (t) in b(0 is defined by, 


e(0 =__ 

b(0 T V noise V„oise T b(0 


( 3 ) 


This is an extension of MUSIC preserving the 
temporal structure. If the contribution of a (t) is 


large in b(^) , then the contribution score e{t) 


becomes large because V noise is almost orthogonal to 

b(0- 


3 Results 

3.1 Experiment 

MEG data of binaural AEF (Auditory Evoked Field) 
related to sound localization brought by ITD 
(Intreaural Time Difference) were analyzed by this 
methods. The stimuli for both ears were computer¬ 
generated pseudo noise (PN). Three types of stimuli 
were designed; 

1. Type A: Silent to correlated pairs of PNs. 

2. Type B: Correlated pairs of PNs to uncorrelated 
pairs of PNs. 

3. Type C: Correlated pairs of PNs to correlated 
pairs of PNs with ITD=135|as (corresponding to 
35 degree). 

Only type A had monaural event factor among 
stimulus. The AEF was measured with a whole-head 



64ch biomagnetometer system (CTF Systems Inc.). 
The sampling frequency was 625Hz. The data were 
100 times averaged and bandpass filtered (2-40Hz). 
Stack plots and iso-contour maps corresponding to 
the stimuli are shown in Fig. 1. 
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Figure 1: MEG related to the binaural stimuli. Each 
column corresponds to the stimuli A-C. The 
subscribed numbers of the maps (in ms) indicate the 
latencies. 

3.2 Spatial characterization 

By performing PCA to the AEF data, the curves of 
the cumulative contribution ratio were obtained (Fig. 
2 ). 



Figure 2: Spatial characterization using PCA. 
Curves of the cumulative contribution ratio are 
graphically presented for the types A-C. Dashed line 
indicates 95%. 

3.3 Temporal characterization 

The contribution of the activity in the data of type A 
in that of type B was calculated using the virtual 
beamformer as shown in Fig. 3. 64 samples in the 
data of type A having AEF activity were selected as 
the activity data (as displayed the rectangle of the 



top column). The contribution score £{t) of the 
activity data was calculated for the data of type B. 
The plot of the contribution score e{t) was obtained 
(as the bottom column). Several peaks of e(t) were 
seen during 100-300ms in latency. 



Figure 3: Temporal characterization using the 
virtual beamformer analysis. Top: Stack plot of 
activity of Type A. Middle: Stack plot of Type B 
(target). Bottom: plot of the contribution score £(t ) . 

4 Discussion 

A typical AEF characterized by activities in both 
temporal lobes can be seen in Fig. 1. In Fig. 2, the 
cumulative contribution ratio reached 94% up to 4th 
component and that of type A is lower than those of 
types B and C from the 1st to 4th component. The 
curve of the cumulative contribution ratio of type A 
is increasing more slowly than types B and C along 
the component axis. This indicates that the 
distribution of the components in type A is more 
spread than types B and C. The number of the event 
factors in type A might be lager than types B and C. 
At least, only type A includes the monoaural event 
factor. 

In Fig. 3, several peaks of £{t ) are clearly seen. 
This indicates that the spatial and temporal activities 
in type B contain those activities expressed by the 
multiple waveforms in the rectangle area of type A. 
Some of event factors in type A might also exist in 
type B. 

The proposed methods will be beneficial to analysis 
for multi-modal activities in the brain. 
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